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重要 的 影响 。 当 前 研究 


大 多 聚焦 于 NDVI 与 气候 因子 的 相关 性 分 析 , 少 数 研究 预测 NDVI 时 忽视 了 考虑 滞后 性 以 及 增加 预报 因子 对 提高 模 
型 精确 度 的 影响 。 基 于 此 ,本 文 比 较 多 元 线性 回归 (MLR)、 人 工 神经 网 络 (ANN ) 和 支持 向 量 机 模型 (SVM ) ,优选 出 


精度 较 高 的 模型 。 在 常规 因子 (降雨 和 气温 ) 基 础 上 增加 了 对 植被 生长 有 影响 的 土壤 湿度 和 日 照 因 子 , 并 考虑 了 不 
同 因子 与 NDVI 之 间 的 滞 时 效应 。 研 究 结果 表明 :(1) 支持 向 量 机 模型 (SVM) 的 拟 合 能 力 最 强 ,NDVI 预 测 精度 最 
高 , 泾 河 . 北 洛 河流 域 均 方 根 误差 均 减 少 1.8% 以 上 。(2) 加 入 土壤 湿度 日照 等 因子 , 泾 河流 域 模型 预测 精度 提高 ， 
泾 河流 域 均 方 根 误差 减少 8.8%。(3) 考虑 清 时 后 , 泾 河北 洛 河流 域 均 方 根 误差 分 别 减少 15% 11% ,NDVI 预 测 精度 
进一步 提高 ,增加 了 模型 的 可 靠 性 。 预 测 结果 可 对 未 来 制定 生态 保护 策略 ,指导 生态 修复 具有 重要 参考 意义 。 
关键 词 : 支持 向 量 机 ， 预测 精度 ; 预测 模型 ， 滞 时 效应 ;植被 覆盖 度 


植被 是 大 气 .土壤 流域 之 间 的 天 然 纽 带 ,是 陆 
地 生态 系统 的 重要 组 成 部 分 ”, 亦 是 气候 圈 、 水 圈 、 
岩石 圈 、 生 物 圈 及 人 类 活动 相互 影响 相互 作用 的 统 
一 体现 。 植 被 覆盖 度 作为 衡量 地 表 植被 状况 的 一 
个 最 重要 的 指标 ,其 在 地 表 植 物 蒸腾 作用 、 光 合作 
用 土壤 水 分 蒸发 等 方面 发 挥 着 重要 作用 。 其 变 
化 作为 土地 覆盖 变化 的 主要 表现 形式 之 一 ,间接 反 
映 了 人 类 活动 对 自然 生态 的 影响 程度 ,这 一 影响 对 
水 循环 过 程 极其 显著 。 归 一 化 植被 指数 (NDVI) 被 
作为 重要 指标 之 一 ,能 有 效 的 表征 地 表 植 被 覆盖 度 
及 植物 生长 趋势 ,其 作为 遥感 近 红 外 波段 与 红外 波 
段 反 射 率 之 差 与 两 者 之 和 的 比值 参数 ,对 植被 覆盖 
状况 的 表现 具有 敏感 性 的 量 测 ,也 被 公认 为 是 表现 
植被 生长 状况 的 有 效 参 数 指标 。NDVI 在 不 同 地 
区 的 差异 化 表现 可 以 反映 这 些 地 区 生态 系统 不 同 
程度 的 生产 力 发 展 状 况 ,不 同年 的 NDVI 可 反映 由 于 
气象 要 素 变化 所 造成 的 系统 生产 力 的 年 际 变化 = 。 
因此 ,NDVI 已 经 被 广泛 应 用 于 植被 时 空 分 布 变 化 
及 演变 的 研究 中 。 植 被 覆盖 指数 在 气候 变化 环境 
下 ,受到 多 种 因素 的 影响 ,呈现 出 复杂 的 非 线性 关 
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系 , 这 增加 了 对 NDVI 准 确 性 预测 的 难度 。 因 此 , 融 
合 多 种 因素 预报 NDVI, 对 提高 NDVI 预 测 的 精确 度 
十 分 重要 。 

目前 ,很 多 学 者 为 了 探究 归 一 化 植被 指数 与 气 
候 要 素 之 间 的 关联 ,对 二 者 进行 了 大 量 的 相关 性 分 
析 研 究 。 杜 加 强 等 “在 一 致 性 检验 基础 上 ,融合 
GIMMS NDVI 数 据 与 MODIS NDVI 数 据 , 将 短 时 间 
序列 拓展 到 长 时 间 序 列 , 分 析 了 气候 变化 对 植被 的 
影响 ,结果 表明 气温 主要 在 春秋 季 促 进 植被 生长 ， 
夏季 降水 量 是 调节 植被 生长 的 主要 因子 。Wu 等 " 
研究 发 现 , 温 度 影响 高 纬度 地 区 植被 覆盖 度 的 大 
小 ,是 植被 生长 的 主要 因素 之 一 。 陶 帅 等 汪 利 用 趋 
势 分 析 和 指数 分 析 等 方法 ,分 析 了 植被 覆盖 度 在 年 
尺度 和 季 尺 度 上 的 时 间 变 化 及 未 来 变化 趋势 ,并 对 
NDVI 时 空 变化 特征 及 相关 因子 进行 交互 作用 分 
析 ,表明 植被 覆盖 度 受 海拔 .年 均 温 .土壤 类 型 及 人 
口 密度 多 种 因子 共同 影响 。 李 登科 等 ”采用 相关 系 
数 .变化 斜率 等 方法 分 析 了 陕 北 长 城 沿线 风沙 区 植 
被 变化 与 降水 .气温 等 气候 因子 之 间 的 关系 ,发 现 
不 同 季 节 NDVI 变 化 幅度 不 一 ,降水 量 是 引起 NDVI 
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年 际 波动 的 主要 因子 ,人 类 活动 等 非 气 候 因 素 也 是 
引起 植被 变化 的 原因 。 对 于 干旱 半 干 量 地 区 的 植被 ， 
NDVI 还 与 土壤 的 特性 ,植被 的 生根 作用 等 有 关 "", 这 
些 因 素 共 同 决定 了 植物 的 生长 状况 及 其 覆盖 度 的 
变化 。 王 伦 澈 5 采用 碳 通 量 数据 对 不 同类 型 植被 
生产 力 模 型 进行 标定 ,与 MOSDIS 算 法 比较 ,分 析 了 
影响 不 同类 型 植被 生长 的 主要 因子 ,发 现 总 辐射 和 
降水 是 影响 研究 区 域 植被 生长 季节 变化 的 关键 因 
子 。 虽 然 以 上 研究 均 表 明 植 被 覆盖 度 受 多 种 因子 
的 共同 作用 ,但 目前 关于 NDVI 预 测 研究 中 多 以 降 
雨 和 气温 作为 输入 因子 。 张 满 固 等 中 基于 降水 和 
气温 数据 ,结合 了 3 种 优化 算法 (网 络 搜 索 法 、 遗 传 
算法 ,粒子 群 算法 ) 对 以 支持 向 量 机 构建 的 NDVI 短 
期 预测 模型 参数 进行 了 优选 ,并 利用 优选 的 最 佳 算 
法 从 不 同 角 度 构建 了 单项 预测 模型 且 进 行 线 型 组 
合 对 NDVI 进 行 回 归 预 测 ,结果 表明 组 合 模型 可 对 
NDVI 进 行 有 效 预测 。Zhou 等 .将 NDVI 作为 因 变 


化 植被 指数 预测 模型 ,这 对 未 来 制定 保护 生态 策 
略 , 指 导 当 地 生态 修复 具有 重要 参考 意义 。 

本 文 在 对 比分 析 了 以 降水 和 气温 数据 构建 的 
支持 向 量 机 (SVM) ”、 人 工 神经 网 络 (ANN) ”和 多 
元 线性 回归 模型 (MLR) ”模拟 植被 变化 的 基础 上 ， 
融合 土壤 湿度 .日 照 等 因子 的 信息 ,构建 了 考虑 不 
同 因 子 与 植被 生长 状况 之 间 沸 后 性 的 NDVI 预测 模 
型 ,以 期 较 好 地 提高 NDVI 预 测 精度 ,为 该 地 区 的 农 
业 生 产 , 旱 灾 预 警 ,生态 发 展 治理 等 方面 提供 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

泾 河 与 北 洛 河 是 渭河 流域 的 两 条 主要 支流 (图 
1) ,处 于 我 国 西北 地 区 的 东部 ,属于 大 陆 性 季风 气 
候 ,降水 主要 集中 在 夏季 ,冬季 两 量 偏 少 ”。 流 域 
内 气温 季节 性 变化 明显 。 流 域 土壤 类 型 主要 为 黄 
土 ,其 土壤 含水 量 较 少 。 受 光照 时 间 较 短 ,年 平均 


tit ,降水 与 蒸 散 发 作为 自 变 量 ,输入 多 元 线性 回归 
模型 预测 NDVI, 结果 表明 基于 降水 和 蒸 散发 的 多 
元 回归 模型 可 以 有 效 预测 不 同 季节 植被 的 动态 变 
化 。 然 而 ,植被 生长 状况 的 优良 程度 并 不 是 只 受到 
单一 气候 因子 的 影响 , 除 降雨 和 气温 ,还 和 多 种 气 
候 要 素 的 变化 有 关 , 受 到 多 种 因子 的 共同 作用 。 基 
于 此 ,有 必要 增加 与 NDVI 密 切 相关 的 更 多 的 输入 
因子 ,这 对 预测 归 一 化 植被 指数 ,提高 模型 可 靠 性 
有 重要 意义 。 

土壤 湿度 作为 重要 的 陆地 生态 系统 输入 因子 ， 
是 水 分 和 能 量 交换 的 重要 因素 ,其 对 植物 作物 的 生 
长 起 到 决定 性 的 作用 中 。Tian 等 中 研究 发 现 , 土 壤 
中 的 含水 量 影响 植被 可 获得 的 储 水 量 ,在 高 纬度 地 
区 ,温度 及 日 照 对 植被 生长 状况 的 驱动 力 比 降水 要 
强 。 裴 志 林 等 “基于 探测 器 模型 对 年 最 大 植被 覆 
盖 度 空间 分 布 的 影响 因子 进行 定量 分 析 ,表明 影响 
植被 覆盖 度 的 主导 因子 是 降水 与 土壤 类 型 和 土地 
利用 方式 的 交互 作用 。 刘 静 用 趋势 分 析 四 检验 、 
皮尔 逊 相关 分 析 与 多 元 线性 回归 模型 从 年 尺度 和 
季节 尺度 分 析 了 不 同 植被 覆盖 区 域 气温 .降水 、 湿 
度 .日照 等 因素 对 NDVI 的 作用 机 制 ,发现 湿度 和 日 
照 对 植被 生长 有 影响 。 基 于 此 ,可 考虑 选取 模型 预 
测 因子 时 增加 土壤 湿度 和 日 照 数据 , 同 时 ,由 于 植 
被 生长 状况 和 和 气候 要 素 之 间 并 非 简单 的 线性 关系 ， 
因而 构建 考虑 涌 时 的 可 靠 的 融合 多 种 因素 的 归 一 


日 照 时 长 约 5h。 主 要 植被 为 落叶 松 ` 山 杏 、 山 桃 、 刺 
槐 及 草地 和 农作物 等 。 将 渭河 流域 的 两 条 主要 文 
流 泾 河和 北 洛 河 作为 研究 区 。 其 中 , 泾 河流 域 灌区 
面积 占 泾 河流 域 面积 的 0.38%。 北 洛 河 农业 灌溉 面 
积 更 大 , 占 北 洛 河流 域 面 积 的 0.66%。 流 域 地 理 位 
置 如 图 1 所 示 。 
1.2 数据 来 源 

采用 的 数据 包括 气象 数据 .土壤 数据 和 NDVI 
数据 ,分 别 来 源 于 中 国 气 象 数据 网 (http:/data.cma. 
cn/) 全 球 陆 面 数据 同化 系统 的 CGLDAS 产品 (http:/ 
disc.sci.gsfc.nasa.gov/hydrology/data-holdings) 及 地 理 
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图 1 泾 河和 北 洛 河流 域 位 置 及 气象 站 点 空间 分 布 图 


Fig. 1 Location and the spatial distribution of meteorological 


stations of Jing River Basin and Beiluo River Basin 
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空间 数据 云 (http://www.gscloud.cn)。 其 中 ,气象 数 
据 为 泾 河流 域 环 县 、 固 原 、 平 凉 、 西 峰 镇 长 武 和 北 
治 河流 域 吴 旗 、 延 安 洛 川 . 铜 川 共 9 个 气象 站 的 
1982 年 1 月 至 2013 年 12 月 共 31 年 的 逐 月 降水 ^ 
温 \ 日 照 等 数据 ,并 利用 泰 森 多 边 形 进行 面 数据 转 
化 。 气 象 站 点 空间 分 布 如 图 1 所 示 。 
1.3 研究 方法 
1.3.1 多 元 线性 回归 ”多 元 线性 回归 是 用 两 个 或 多 
个 自 变量 解释 因 变 量 的 一 种 线性 计算 模型 。 多 元 
线性 回归 模型 假设 因 变 量 与 自 变 量 间 存 在 线性 关 
系 ,用 最 小 二 乘法 ,计算 得 到 每 个 自 变量 对 应 的 
回归 系数 ,建立 多 元 线性 方程 ,对 自 变 量 实现 预 
报 。 一 般 多 元 线性 回归 方程 如 下 : 

Y-B,- p,X,* B,X,*--- B, X, te (1) 
式 中 :7 为 因 变 量 ; B, 为 自 变 量 对 应 的 回归 系数 ,n= 
0,1,…,m; 训 是 自 变量 ,n=0,1,…,m; e 则 是 消除 多 
个 自 变 量 对 因 变 量 影响 的 随机 误差 值 。 
1.3.2 支持 向 量 机 文 持 向 量 机 是 通过 确定 核 函 
数 ,构建 最 优 超 平面 ,解决 了 复杂 的 非 线性 拟 合 问 
题 。 对 已 知 的 数据 集 , 文 持 向 量 机 通过 非 线性 映 
射 ,将 数据 映射 到 高 维特 征 空间 ,利用 线性 问题 求 
解 方法 在 高 维特 征 空间 寻求 最 优 的 回归 函数 3 。 
其 中 ,线性 回归 函数 如 公式 (2): 

f(x)=w, x)+b (2) 

WP: fe) 为 非 线性 映射 函数 ; w 为 权重 ;x 为 第 i 个 
样本 的 输入 变量 九 为 闽 值 。 
1.3.3 人 工 神 经 网 络 人 工 神经 网 络 是 一 种 反 向 误 
差 传 播 ,与 动物 神经 网 类 似 , 不 断 修 正 更 新 权 值 和 
浆 值 的 方法 ,用 误差 反 向 传播 调整 网 络 参 数 ” ,是 
一 种 被 广泛 应 用 具有 较 强 的 非 线性 映射 能 力 和 自 
学 习 能 力 的 算法 后 。 计 算 公 式 如 下 : 


n-l 
xe Sos-0) (70, 1, 2, =, n=) G) 
k=0 


| (n- 0, 1; 2, kag. m- 1) (4) 
i=0 


AP Uc TIO HS Aa A SRE s wa RIA 
入 层 第 大 个 神经 元 到 隐 含 层 第 ;个 神经 元 权 值 ; ww 
表示 输入 值 ;w 表 示 从 隐 含 层 第 :个 神经 元 到 输出 
层 第 n 个 神经 元 的 权 值 ; 9 为 隐 含 层 第 ;个 神经 元 国 
值 ; 0, 为 输出 层 第 "个 神经 元 阔 值 。 

人 工 神 经 网 络 由 输入 变量 、 输 出 变量 、 隐 含 层 
及 激励 函数 组 成 。 学 习 过 程 中 ,首先 对 样本 进行 预 
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处 理 , 确 定 网 络 结构 , 设 定 权 值 及 阅 值 。 归 一 化 处 
理 数 据 后 ,对 网 络 进行 训练 ,模拟 网 络 ,得 到 权 值 阔 
值 的 最 优 值 ,计算 实际 输出 值 及 目标 输出 间 的 误 
2 ,检验 误差 ,是否 达 到 预 设 精度 ,达到 则 输出 结 
A ,未 达到 则 重新 模拟 网 络 ,计算 误差 。 人 工 神经 
网 络 具 有 较 强 针对 非 线 性 映射 的 能 力 , 有 较 好 的 自 
适应 学 习 能 力 和 较 强 的 鲁 棒 容 错 等 特性 ” 。 

1.3.4 考虑 滞 时 的 多 因子 预测 模型 ”为 了 使 预测 模 
型 更 可 靠 ,提高 NDVI 预 测 精 度 ,本 文 考 虑 了 因子 的 
滞后 性 ,并 通过 添加 除 降 雨 .气温 以 外 的 其 他 因子 
(土壤 湿度 和 日 照 ) 的 方式 ,构建 了 多 因子 预测 模 
型 ,以 期 提高 预测 模型 的 精度 。 

首先 , 仅 将 月 降雨 和 月 气温 数据 作为 输入 变 
ti ,分别 输入 至 BP 神 经 网 络 .支持 向 量 机 与 多 元 线 
性 回归 模型 中 比较 其 预测 性 能 ,并 优选 出 最 佳 预测 
模型 。 其 次 ,增加 土壤 湿度 和 日 照 月 数据 作为 输入 
因子 ,与 只 输入 降雨 及 气温 的 预测 结果 进行 对 比 ,分 
析 预 测 精 度 是 否 有 所 提高 。 最 后 ,融合 NDVI 与 降雨 、 
气温 .土壤 湿度 .日照 因子 的 清 时 数据 并 进行 预测 。 

本 文选 用 纳什 系数 (NSE) 2” 、 均 方 根 误差 
(RMSE) 所 和 平均 相对 误差 (MRE) 2 构建 评价 指标 
体系 ,评估 预测 精度 。 其 中 ,纳什 系数 验证 模型 模 
拟 结果 的 好 坏 ,计算 公式 如 下 : 

Yi 
NSE=1- 人 (5) 
2.6; -7) 
IUP y KA WE s c; Js ME; 了 是 观测 值 的 平均 
值 。 纳 什 系数 趋 近 于 1, 表示 模型 结果 好 ,可 信和 度 
高 , 趋 于 0, 则 模型 结果 差 , 可 信和 度 低 。 

均 方 根 误差 表示 观测 值 与 真实 值 之 间 误 差 的 

指标 ,计算 公式 如 下 : 


y 入 一 二 
二 -一 一 (6) 


n 
式 中 :y 是 观测 值 ;# 是 预测 值 ,n 为 样本 数 。 均 方 根 
误差 值 越 小 ,真实 值 越 接近 观测 值 。 

平均 相对 误差 能 很 好 的 反应 预测 值 误差 的 实 
际 情况 ,计算 公式 如 下 : 

MRE = 二 (7) 

式 中 :y 是 观测 值 ;# 是 预测 值 ,n 为 样本 数 。 流 程 图 
如 图 2 所 示 。 


RMSE = 
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预测 结果 分 析 


图 2 流程 图 
Fig.2 Flow diagram 


2 结果 与 分 析 


2.1 确定 最 优 预测 模型 

将 1982—2013 年 泾 河流 域 月 降雨 和 气温 数据 
作为 输入 向 量 , 分 别 用 多 元 线性 回归 模型 人工 神 
经 网 络 模 型 及 支持 向 量 机 模型 对 研究 区 归 一 化 植 
被 指数 (NDVI) 进 行 预测 ,训练 期 为 1982 年 1 月 至 
2006 年 12 月 ,验证 期 为 2007 年 1 月 至 2013 年 12 月 ， 
输出 预测 结果 ,不 同 模型 对 NDVI 的 预测 精度 各 不 
相同 。 比 较 纳什 效率 系数 、 均 方 根 误差 及 平均 相对 
误差 ,对 于 不 同 模型 的 NDVI 预 测 结果 如 表 1 所 示 。 

在 泾 河流 域 , 用 多 元 线性 回归 模型 预测 的 结 
AR ,验证 期 其 纳什 系数 仅 为 0.731, 均 方 根 误差 为 
0.092 ,平均 相对 误差 为 0.224; 采 用 人 工 神 经 网 络 模 
型 计算 的 结果 ,训练 期 纳什 系数 比 多 元 线性 回归 模 
型 高 2.53%; 而 用 支持 向 量 机 模型 预测 的 结果 ,验证 
期 纳什 系数 最 大 ,为 0.761 ,训练 期 纳什 系数 比 人 工 
神经 网 络 高 1.47% , 比 多 元 线性 回归 模型 高 4.04% 。 
因此 ,在 泾 河流 域 ,支持 向 量 机 模型 比 多 元 线性 回 
归 模 型 和 人 工 神 经 网 络 模型 更 优 , 预 测 精度 更 高 。 


在 北 洛 河流 域 , 多 元 线性 回归 模型 预测 的 结果 , 验 
证 期 纳什 系数 为 0.544, 均 方 根 误差 0.128, 用 人 工 神 
经 网 络 模 型 计算 ,验证 期 纳什 系数 为 0.611, 均 方 根 
误差 0.120。 支 持 向 量 机 模型 预测 结果 比 多 元 线性 
回归 模型 预测 结果 高 21.9% , 比 人 工 神 经 网 络 模型 
预测 结果 高 8.52% 。 训 练 期 支持 向 量 机 模型 预测 
NDVI 纳 什 系数 与 多 元 线性 模型 相 比 ,高 14.64% , 支 
持 向 量 机 模型 预测 结果 纳什 系数 比 人 工 神经 网 络 
模型 高 0.01%。 总 体 上 支持 向 量 机 模型 仿真 模拟 结 
果 最 好 。 

结果 表明 ,3 种 模型 中 ,在 处 理 复杂 的 非 线性 映 
射 关 系 时 ,人 工 神经 网 络 模型 可 利用 已 有 的 数据 资 
料及 经 验 , 通 过 复杂 的 表达 方式 ,表现 出 输入 变量 
及 输出 变量 之 间 的 非 线 性 关系 。 人 工 神经 网 络 模 
型 比 多 元 线性 回归 模型 的 预测 能 力 更 好 ,有 更 优 的 
泛 化 能 力 。 而 在 同一 环境 下 ,输入 相同 因子 时 , 支 
持 向 量 机 的 优势 在 于 把 优化 问题 转化 为 凸 二 次 规 
划 问 题 , 而 人 工 神 经 网 络 模型 在 得 全 局 最 优 解 时 会 
出 现 局 部 极 值 的 问题 ,支持 向 量 机 模型 避免 了 这 一 
缺点 , 且 相 较 于 多 元 线性 回归 模型 ,其 还 可 以 解决 
复杂 的 非 线 性 问题 。 故 支持 向 量 机 模型 的 预测 精 
度 比 前 两 者 更 高 。 
2.2 预测 模型 输入 数据 的 选择 

大 多 研究 只 考虑 降雨 或 气温 单 因子 对 NDVI 的 
影响 ,但 除 降雨 和 气温 ,NDVI 还 和 多 种 气候 要 素 的 
变化 有 关 , 受 到 多 种 因子 的 共同 作用 。 如 图 3 所 示 ， 
NDVI 与 不 同 气 候 因 子 和 土壤 湿度 之 间 均 存在 复杂 
的 非 线 性 关系 。 故 为 提高 模型 模拟 预测 结果 的 准 
确 性 ,有 必要 增加 与 NDVI 密 切 相 关 的 更 多 的 输入 
因子 。 本 研究 将 土壤 湿度 和 日 照 加 入 降雨 气温 模 
型 的 输入 变量 中 ,预测 NDVI, 将 其 预测 结果 与 降雨 
气温 作为 输入 变量 的 同一 模型 的 预测 结果 进行 比 
较 , 由 前 述 研 究 发 现 支持 向 量 机 模型 预测 精度 更 
高 ,故此 选用 支持 向 量 机 作为 预测 模型 。 预 测 结 


表 1 多 元 线性 回归 、 人 工 神 经 网 络 及 支持 向 量 机 对 NDVI 预 测 


Tab.1 Multiple linear regression, artificial neural network and support vector machine predict NDVI 


多 元 线性 回归 
人 研究 区 域 y 
" 训练 期 0.820 0.066 0.203 0.936 
泾 河 u 
验证 期 ” 0.731 0.092 0.224 0.750 
o 训练 期 0.832 0.077 0.187 0.953 
JRI 0 
验证 期 0.544 0.128 0.205 0.611 


SERES 支持 向 量 机 


别 Ã—— n ———M—————————— ——M——————— 
纳什 系数 均 方 根 误差 平均 相对 误差 ”纳什 系数 均 方 根 误差 平均 相对 误差 ”纳什 系数 均 方 根 误差 平均 相对 误差 


0.038 0.118 0.928 0.040 0.115 
0.088 0.183 0.761 0.086 0.177 
0.040 0.092 0.953 0.040 0.083 
0.120 0.191 0.663 0.103 0.184 
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图 3 NDVI 与 各 因子 相关 性 分 析 结 : 
Fig.3 Results of correlation analysis between NDVI and various factors 
如 表 2 所 示 。 型 精度 与 多 因子 预测 模型 精度 相差 较 小 。 


由 表 2 可 知 ,在 泾 河 ,输入 变量 增加 土壤 湿度 和 
日 照 因 子 后 ,预测 精度 提高 。 只 输入 降雨 和 气温 两 
因子 预测 模型 的 纳什 系数 为 0.761 ,而 多 因子 输入 变 
量 模型 预测 结果 ,其 均 方 根 误差 比 两 因子 预测 模型 
减少 了 8.806% ,平均 相对 误差 减少 5.649% , 总体 上 
多 因子 预测 模型 精度 比 两 因子 预测 模型 精度 要 
高 。 在 北 洛 河流 域 , 降 雨 气温 两 因子 预测 模型 纳什 
系数 为 0.663 ,多 因子 预测 模型 纳什 系数 比 两 因子 预 
测 模型 纳什 系数 高 1.06% ,降雨 气温 两 因子 预测 模 


2.3 区 域 植被 对 气候 因子 的 滞后 性 

Wu 等 "利用 神经 网 络 算法 ,构建 了 NDVI- 降 雨 
模型 ,其 结果 验证 了 NDVI 与 降雨 之 间 存 在 滞后 性 ， 
与 生态 规律 一 致 。 故 植被 生长 状况 并 非 只 取决 于 
现在 的 气候 条 件 ,与 前 期 气候 条 件 状 况 有 很 大 的 关 
系 , 所 以 本 文 计算 了 各 个 气候 因子 .土壤 湿度 与 ND- 
VI 之 间 的 相关 系数 。 

由 不 同 区 域 各 因子 与 NDVI 之 间 的 相关 系数 可 
知 ( 表 3) ,植被 变化 对 气温 .降水 .土壤 湿度 及 日 照 


表 2 多 因子 输入 模型 预测 结果 比较 


Tab. 2 Comparison of prediction results in multiple factors input-model 


MIS f 降雨 气温 两 因子 降雨 气温 .土壤 湿度 日照 多 因子 
研究 区 域 模型 类 别 - — — — - — — — 
纳什 系数 均 方 根 误差 平均 相对 误差 纳什 系数 均 方 根 误差 平均 相对 误差 
-— 训练 期 0.928 0.040 0.115 0.942 0.036 0.099 
m 
E 验证 期 0.761 0.086 0.177 0.729 0.079 0.167 
u 训练 期 0.953 0.040 0.083 0.959 0.041 0.084 
北 洛 河 
验证 期 0.663 0.103 0.184 0.670 0.118 0.188 
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33 不 同 研究 区 域 各 气候 因子 与 NDVI 之 间 相 关系 数 


Tab.3 Correlation coefficient between various climatic factors and NDVI in different regions 


研究 区 域 气候 因子 "en 
0 1 2 3 4 5 6 7 8 9 10 11 12 
泾 河 NDVI - 0.80 0.42 0.01 -0.42 -0.78 -0.90 -0.78 -0.41 0.01 0.42 0.78 0.97 
降雨 0.90 0.88 0.62 0.19 -0.30 -0.70 -0.91 -0.88 -0.62 -0.18 0.31 0.69 0.90 
气温 0.68 | 0.60 032 -0.04 -0.35 -0.59 -0.69 -0.61 -037 0.01 0.37 0.59 0.68 


土壤 湿度 -0.15 -0.32  -041  -038 -0.23  -0.05 0.15 0.33 0.43 0.41 0.26 0.03 -0.20 
g 0.20 0.33 0.35 0.31 0.17 0.01 -0.17 -0.30 -0.38 -0.36 -0.24  -0.02 0.18 


北 治 河 NDVI 一 0.42 0.00  -043 -0.79 -0.92 -0.78 -0.42 0.01 0.43 0.79 0.98 0.79 
降雨 0.90 | 087 0.61 0.17  -031 -0.71 -0.92 -0.88 -0.61 -0.16 032 071 0.90 
气温 0.75 | 0.64 033 -0.08 -0.41 -0.64 -0.73 -0.63 -0.38 003 042 0.66 075 
土壤 湿度 -0.20  -037 -0.45 -0.40 -0.22 -0.01 0.17 0.33 0.42 0.40 0.26 0.001 | -022 
H RH 0.6 |— 030 0.32 029 015 -0.01 -0.15 -0.25 -0.32 -0.31 -021 -0.02 014 
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等 因子 的 滞后 时 间 是 不 同 的 ,这 与 一 些 学 者 的 研究 
结果 相符 ,Anderson 等 研究 了 不 同 气候 因子 的 滞 
后 性 ,发 现 植被 对 气候 因子 存在 响应 时 间 , 且 对 气 
候 因 子 响 应 的 滞 时 不 同 。 因 此 ,在 探讨 气候 一 植被 
交互 作用 机 制 时 ,应 考虑 灌 后 性 。 选 用 气候 因子 与 
NDVI 之 间 绝 对 值 大 于 0.3 的 相关 系数 对 应 的 月 作 
为 预选 滞 时 ,相关 系数 满足 95% 置 信和 度 检 验 。 
此 , 泾 河流 域 NDVI 滞 时 是 1.2.4.5.6.7、.8、10、11、 
124 H ; ERE WERTE E 0.1,2,5,6,7,8,10,11,12 
月 ;土壤 湿度 滞 时 为 1.2.3. 7.8 .9 个 月 ;日 照 滞 时 选 
为 1.2.3.8.9 个 月 。 同 理 ,优选 得 到 北 洛 河流 域 
NDVI 与 各 气候 因子 的 滞 时 。 

将 上 述 优选 结果 考虑 滞 时 的 NDVI 降雨 、 气 
温 .土壤 湿度 .日 照 各 因子 作为 输入 项 ,以 输入 因 
子 .目标 因子 建立 支持 向 量 机 模型 ,将 其 预测 结 
与 未 考虑 清 后 时 间 影 响 的 预测 结果 进行 比较 ( 表 
4)。 泾 河流 域 , 考 虑 滞 时 的 预测 模型 其 结果 ,模型 
的 均 方 根 误差 比 未 考虑 灌 时 模型 的 减 小 24.39% , 平 
均 相 对 误差 减 小 23.50% ;在 北 洛 河流 域 ,考虑 浪 后 
时 间 的 模型 预测 结果 纳什 系数 比 不 考虑 浪 后 时 间 


的 模型 结果 高 34.93%, 均 方 根 误差 减 小 51.3%, 平 均 
相对 误差 减 小 39.46%。 对 泾 河流 域 和 北 洛 河流 域 
训练 期 .验证 期 考虑 滞 时 模型 与 不 考虑 清 时 模型 的 
预测 结果 进行 比较 (图 4) ,考虑 各 因子 .土壤 湿度 与 
NDVI 之 间 潍 时 的 模型 预测 精度 比 未 考虑 汪 时 的 模 
型 预测 精度 显著 提高 。 因 此 ,变化 环境 下 考虑 滞后 
时 间 的 模型 对 植被 生长 状况 预测 是 有 潜力 的 。 


3 讨论 


上 述 结果 表明 ,考虑 滞 时 的 预测 模型 精度 更 
高 。 这 可 能 是 由 于 植物 本 身 有 一 定 的 耐 受 力 ,前 期 
内 部 有 定量 的 水 分 储存 ,在 干旱 发 生 时 ,植被 有 相 
应 的 应 对 策略 ,自身 应 急 机 制 啊 应 ,对 气候 变化 的 
发 生存 在 反应 时 间 , 故 其 不 能 对 气候 变化 立即 做 出 
回应 ,其 生长 状况 受到 前 期 气候 要 素 变化 的 影响 较 
大 。 本 研究 在 以 往 影 响 植被 生长 的 常规 因子 的 研 
究 中 增加 了 土壤 湿度 .日 照 等 敏感 因子 ,结果 表明 ， 
在 泾 河流 域 降雨 .气温 基础 上 融合 了 土壤 湿度 和 日 
照 数 据 后 模型 预测 精度 提高 ,该 结果 可 能 与 植被 生 


表 4 支持 向 量 机 模型 考虑 /不 考虑 灌 时 的 预测 结果 


Tab.4 Prediction results based on support vector machine with/without time lag 


考虑 滞 时 的 预测 模型 


不 考虑 浪 时 的 预测 模型 


研究 区 域 类 别 


纳什 系数 均 方 根 误差 平均 相对 误差 纳什 系数 均 方 根 误差 平均 相对 误差 
— 训练 期 0.982 0.020 0.048 0.940 0.037 0.103 
ET. 
T 验证 期 0.880 0.062 0.140 0.787 0.082 0.183 
a 训练 期 0.989 0.020 0.038 0.947 0.042 0.085 
北 洛 河 - 
验证 期 0.927 0.056 0.112 0.687 0.115 0.185 
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图 4 泾 河北 洛 河 流域 训练 期 和 验证 期 NDVI 的 模拟 结 


Fig. 4 Simulation results of NDVI in the training period and verification period of Jing River Basin and Beiluo River Basin 


长 状况 的 影响 因素 并 不 单一 有 关 : 植 被 覆盖 状况 在 
气候 变化 环境 下 , 受 多 种 因素 的 影响 。 而 在 北 洛 
河 ,预测 精度 并 无 明显 变化 ,可 能 是 由 于 该 研究 区 
域 有 灌溉 面积 , 泾 河 研 究 区 灌区 面积 占 流域 面积 的 
0.38%。 北 洛 河 农业 灌溉 面积 更 大 , 占 流域 面积 的 
0.66%, 受 人 类 活动 影响 较 大 。 由 于 不 同人 研究 区 域 
NDVI 敏感 因子 不 同 , 所 以 本 研究 结果 虽 不 一 定 适 
用 于 所 有 区 域 ,但 本 人 研究 运用 了 多 元 线性 回归 .BP 
神经 网 络 模型 和 SVM 模型 3 种 预测 模型 ,这 3 种 模 
型 对 研究 区 域 的 选择 无 要 求 , 各 个 区 域 均 可 适用 ， 
所 以 不 同 区 域 的 研究 ,需要 优选 敏感 因子 ,但 3 种 预 
测 模型 和 框架 可 推广 到 其 他 区 域 。 通 过 目前 发 生 
的 气候 变化 ,为 准确 预测 下 一 阶段 植被 生长 状况 的 
变化 情况 ,制定 应 对 策略 ,指导 生态 修复 达到 保护 
生态 目的 提供 理论 依据 。 下 一 阶段 可 从 3 个 方面 突 
破 :(1) 对 NDVI 的 影响 机 理 深 入 研究 ,融合 遥 相 关 
因子 ,对 输入 变量 实现 预 处 理 。(2) NDVI 预 测 模型 
的 进一步 优化 ， a A 
高 模型 预测 精度 。(3) 定量 分 解 气候 变化 与 人 类 
动 对 NDVI 变 化 的 影响 。 


4 结 论 


本 文 将 降雨 .气温 作为 输入 变量 ,分 别 采用 多 
元 线性 回归 模型 人 工 神 经 网 络 模型 和 支持 向 量 机 
模型 对 渭河 流域 的 支流 ( 泾 河北 洛 河流 域 ) 的 ND- 
VI 进行 了 预测 ,通过 分 析 3 种 模型 的 性 能 优选 出 适 
合 研 究 区 的 最 佳 预测 模型 ,并 融合 土壤 湿度 和 日 照 
数据 作为 输入 因子 。 同 时 在 分 析 了 各 输入 因子 的 


滞后 性 对 植被 生长 覆盖 程度 的 影响 基础 上 构建 了 
研究 区 NDVI 考 虑 清 时 的 预测 模型 ,得 出 如 下 结论 : 

(1) 3 种 模型 的 预测 精度 由 高 到 低 依次 为 : 文 持 
向 量 机 模型 > 人 工 神经 网 络 > 多 元 线性 回归 模型 , 表 
明和 气象 要 素 与 植被 覆盖 度 之 间 存 在 复杂 的 非 线性 
关系 , 且 占 较 大 比重 。 

(2) 文 持 向 量 机 模型 中 ,输入 因子 加 入 土壤 湿 
度 和 日 照 数据 ,精度 比 降雨 、 气温 两 因子 预测 精 
更 高 ,模型 更 可 靠 。 表 明 植 被 生长 状况 并 非 受 单一 
因素 影响 ,而 受 多 种 因素 共同 作用 。 北 洛 河流 域 由 
于 受 人 类 活动 的 影响 ,精度 较 无 显著 变化 。 

(3) 在 对 泾 河北 洛 河 研究 区 植被 归 一 化 指数 
进行 模拟 预测 中 , 相 较 于 输入 无 清 时 的 因子 ,考虑 
滞 时 的 预测 因子 模型 预测 精度 更 高 。 植 被 受 不 同 

气候 因素 的 影响 具有 不 同 的 滞 时 ,考虑 植被 主要 敏 
感 因 子 的 浏 后 性 对 揭示 植被 生长 对 气候 因素 的 反 
应 极其 重要 。 
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Simulation of vegetation change based on BP-SVM mode 
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Abstract: Vegetation is an important link that connects the biosphere, atmosphere and hydrosphere, and has an 
important impact on the watershed in the ecological environment and on the exchange of water and heat. 
Previous studies focused on correlation analyses between the normalized difference vegetation index (NDVI) and 
climate factors, but the lag and the increase of forecast factors were rarely used to improve the precision of the 
NDVI model. Thus, this article compared the prediction performance of the multiple linear regression, artificial 
neural network and support vector machine, then selected the model with the highest precision. Using 
conventional factors (rainfall and temperature), the prediction performance was tested when soil moisture and 
sunshine duration were increased, which affect vegetation growth, and the time lag effect between the different 
factors and NDVI were considered. (1) SVM had the strongest fitting ability and the highest NDVI prediction 
accuracy. The root-mean-square error of the Jing River Basin and the Beiluo River Basin were reduced by more 
than 1.896. (2) The root- mean- square error (RMSE) of Jing River Basin decreased by up to 8.8% when soil 
moisture, sunshine and other factors were added. (3) Considering the lag time, the root-mean-square error of the 
Jing River Basin and the Beiluo River Basin decreased by 15% and 1196, respectively, and the predicting 
accuracy of NDVI was improved further, thus increasing the reliability of the model. These prediction results 
have an important reference significance for the formulation of ecological protection strategies and the guidance 
of ecological restoration in the future. 


Keywords: support vector machines; prediction accuracy; prediction model; time-lag effect; vegetation coverage 


